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ABSTRACT 

We study the H2 photodissociation regions around OB stars in primordial gas clouds whose 
virial temperatures are between a few hundred and a few thousand Kelvin. In such small 
objects, a single O star can photodissociate a mass equal to that of the cloud itself. As a result, 
the clouds deplete their molecular coolant and cannot cool in a free-fall time, and subsequent 
star formation is totally quenched. This indicates that stars do not form efficiently in small 
objects and that these objects contribute little to the reionization of the universe. 



Subject headings: cosmology: theory — H II region — ISM: clouds — molecular processes 
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1. Introduction 

After standard recombination at z ~ 1100, the universe is considered to be reionized at some redshift 
z :> 5 from the negative results of Gunn-Peterson experiments in quasar spectra. It has been pointed out 
that first stars can play an important role in the reionization of the universe (e.g., Couchman & Rees 1986). 
This scenario has been investigated in detail using numerical simulations (e.g., Gnedin & Ostriker 1997) 
or semi-analytical models (e.g., Fukugita & Kawasaki 1994; Haiman & Loeb 1997). In the latter models, 
it is crucial to know whether a cloud of mass M that virializes at redshift z can cool or not, i.e., if star 
formation occurs or not. 

Haiman, Thoul, & Loeb (1996) investigated this problem (see also Tegmark et al. 1997). They 
found that molecular cooling plays a crucial role for clouds with T vir < 10 4 K ("small" pregalactic clouds, 
hereafter) at 10 < z vil <J 100, where T vir and z v j r are the virial temperature and the redshift at virialization 
respectively, and determined the minimum mass of the virialized cloud must have in order to cool in a 
Hubble time. 

On the other hand, Haiman, Rees, & Loeb (1997) pointed out that in the presence of ultraviolet (UV) 
background radiation at the level needed to ionize the universe, molecular hydrogen is photodissociated by 
far ultraviolet (FUV) photons, whose radiation energy is less than the Lyman limit, in small pregalactic 
clouds. Thus, they asserted that molecular hydrogen in small pregalactic clouds is universally destroyed 
long before the reionization of the universe. In their reionization model, Haiman & Loeb (1997) assumed 
that only objects with virial temperature above 10 K can cool in a Hubble time, owing to atomic cooling 
and star formation occurs subsequently. 

Recently, Ciardi, Ferrara, & Abel (1999) found that the photodissociated regions are not large enough 
to overlap at z ~ 20 — 30. In the same redshift range, the flux of FUV background is well below the 
threshold required by Haiman et al.(1997) to prevent the collapse of the clouds. 

However, molecular hydrogen in a virialized cloud is photodissociated not only by external FUV 
background radiation, but also by FUV photons produced by massive stars within the cloud. Here, we 
assess the negative feedback of massive star formation on molecular hydrogen formation in a primordial 
cloud. 



2. Region of Influence of an OB star 

Around an OB star, hydrogen is photoionized, and an HII region is formed. Ionizing photons hardly 
escape from the HII region, but photons whose radiation energy are below the Lyman limit can get away. 
Such FUV photons photodissociate molecular hydrogen, and a photodissociation region (PDR) is formed 
just outside the HII region. In this section, we study how much mass in a primordial cloud is affected by 
such FUV photons from an OB star and, as a result, becomes unable to cool in a free-fall time owing to the 
lack of the coolant. 

We consider a small pregalactic cloud of primordial composition. In such an object, H2 is formed 
mainly by the H~ process at z <^ 100: 

H + e- -> H"+7; (1) 
H + H- -» H 2 +e~. (2) 

At z ;> 100, H~ is predominantly photodissociated by CMB photons before the reaction (J2J) proceeds. On 
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the other hand, in a PDR, photodissociation of H by UV radiation from an OB star does not dominate 
the reaction (||) except in the vicinity of the star. Then we neglect the photodissociation of H~. 

The rate-determining stage of the H~ process is the reaction ([[]), whose rate coefficient fejj- is (de Jong 
1972) 

&H- = 1.0 x 1(T 18 T s _1 cm 3 . (3) 
In a PDR, H2 is dissociated mainly via the two-step photodissociation process: 

H 2 + 7 -> H* ^ 2H, (4) 

whose rate coefficient fc2step is given by (Kepner, Babul, & Spergel 1997; Draine & Bertoldi 1996) 

fc 2s tc P = 1.13 x 10 s F LW s-\ (5) 

Here Flw (ergs s _1 cm _2 Hz _1 ) is the averaged radiation flux in the Lyman and Werner (LW) bands and 
can be written as 

FlW = ^LW, ex/shield j (6) 

where FLw.cx is the incident flux into the PDR at 12.4 eV and the shielding factor /shield is given by (Draine 
& Bertoldi 1996) 



/shield = min 



1,( : 



(7) 



1Q 14 

The timcscale in which the H2 fraction reaches the equilibrium value is given by 

^dis = k 2s l ap - (8) 

In the presence of FUV radiation, if the temperature and density were fixed, the H 2 fraction / initially 
would increase and reach the equilibrium value for a temporal ionization fraction after ~ i<jis) an d then 
it would decline as ionization fraction decreased as a result of recombination. Actually, if the pregalactic 
cloud can once produce a sufficient amount of molecular hydrogen to cool in a free-fall time, the cloud can 
collapse and star formation occurs subsequently. Note that because of their low ionization degree, inverse 
Compton cooling by CMB photons is not effective in small objects with T V1I <J 10 4 K. 

We investigate here how much mass around an OB star is affected by the photodissociating FUV 
radiation from the star and becomes unable to cool in a free-fall time. In particular, we seek the lower 
bound of such mass. 

The equilibrium number density of H2 under ionization degree x is 

fe H - 2 



7lH 2 = T xn (9) 

n<2step 

= 0.88 x lQ- 26 xF~^Tn 2 . (10) 



Near the star, FLw is so large that the dissociation time idis is smaller than the recombination time 
£rcc = (fcrcc^i^) -1 , where fc roc is the recombination coefficient, xi is the ionization degree at virialization, 
and n is the number density of hydrogen nuclei. Then the chemical equilibrium between above processes is 
reached before significant recombination proceeds. Far distant from the star, the recombination proceeds 
before the molecular fraction reaches the equilibrium value and the ionization degree significantly diminishes. 
However, since we are seeking how much mass is at least affected by the photodissociating FUV radiation 
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from the star, we use the equilibrium value (|^) with initial ionization degree x — x\ as the H2 number 
density. 

Consider an OB star, which radiates at the rate of I/lw [ergs s _1 Hz -1 ] in the LW bands. For an 05 
star, whose mass is 40 Mq, Llw — 10 24 ergs s _1 Hz -1 . At the point whose distance from the star is r, the 
averaged flux in the Lyman and Werner bands is approximately given by 

FhW = /shield- (11) 

Using above relations, we obtain the H2 column density between the star and the point whose distance 
from the star is r: 

N = / CV (N H2 < 10 14 cm- 3 ) 

Ha \ 10 14 [0.25(CT/10 14 ) + 0.75] 4 (A^h 2 > 10 14 cm- 3 ) 

where 



C = 0.88 x 10~ 26 XiL^Tn 2 , (13) 
V = f r 3 . (14) 

Here we have assumed n =const. in space, for simplicity. We define here the region of influence around a 
star as that where the cooling time i coo i = (V 2 ) kT b ecomes larger than the free-fall time tg = ( „^ nh ) 1 ^ 2 
as a result of the photodissociation of molecular hydrogen, where / = nn 2 /n is the H2 concentration, Ar- 2 
is the cooling function of molecular hydrogen, and is the baryon mass fraction. 

The condition t coo \ > tg is satisfied as long as the H2 fraction 

f<f (cool) _ f 24Gm P U/2 kT (.r.) 

where we used our fit to the Martin, Schwarz,& Mandy (1996) H2 cooling function 

A H2 ^ 4 x 10- 25 (— ^) 4 ergs s^cm 3 (17) 

for the low temperature (600K <; T <; 3000K) and low density (n < 10 4 cm" 3 ) regime. The same condition 
leads to the condition on the averaged flux in the LW bands with equation M: 



r (cooi) _ k H -Xjn 

1.13 x 10 8 / (cool) 



^LW > ^ LW - 1 -i o . . i»«Wrnnn 



Corresponding to this critical LW flux F^™ 1 ' , a critical radius A 000 ^ is determined by the relation 

F LW [r( cool )] = ^ ol) . (20) 
Actually, the timescale for molecular hydrogen to reach the equilibrium value at r = r( c ° o1 ' , 

= 4.» * jov^-^r^^r'O-"', < 21 » 
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is longer than the lifetime of an OB star £ob ~ 3 x 10 6 yr. This means that as far region as r — r( c ° o1 ) is 
rarely affected within the lifetime of a single massive star. We define here another critical LW flux F^^ 
and its corresponding radius A ca ^ where the timescale for molecular hydrogen to reach the equilibrium 
value £dis becomes equal to the lifetime of an OB star; 



F$ = 0.93 x lCT 22 ergs s ~ W^z^^-)- «22 , 



F VN [r^] = F^>. (23) 

In the region of influence, we require that two conditions be met: (1) the cooling time for the 
equilibrium H2 fraction is longer than the free-fall time (i.e., t coo \[f — / < - cq - ) ] > iff, where /( cc *) is the 
equilibrium H2 fraction) and (2) the equilibrium H2 fraction is reached within the lifetime of the central 
star (i.e., idis < £ob)- Hence the radius of influence r( mf ) is determined by the smaller one of either r( c ° o1 ) 
or r( cq ) ; 

r< inf > = min[r( co °V (cq) ]- (24) 

Note the LW flux Flw at which tdi S becomes equal to t lec is 2 x 
10 _24 (a;i/10~ 4 )(n/lcm _3 )(T/10 3 K) _0 - 64 ergs s _1 cm _2 Hz _1 . Here we used the recombination coefficient 
k X ec = 1-88 x io _10 T -0 ' 64 s~ 1 cm 3 (Hutchins 1976). Therefore, even at the edge of the region of influence, 
the condition idis < ircc is usually satisfied. In this case, the chemical equilibrium value for H2 is reached 
before the ionization degree significantly decreases from the initial value X[. 

If the self-shielding of LW band photons could be neglected, the radius r* 000 ') would be given by 

(Cool) _ (Cool) _ r L LW l/2 

r - r no-sh - L „(cool) i 

This expression is valid only when ]Vh 2 < 10 14 cm -2 . 

When the H2 column density N^ 2 becomes larger than 10 14 cm -2 , self-shielding of LW band photons 
begins as can be seen from equation (fm. Here, we define the shielding radius r s h as the radius where 
-^H 2 ( r sh) = 10 14 cm~ 2 . Using equation dig), the shielding radius is 

- = ijMm f " <27) 
. »»^^[ ,^, w^rn^. (28) 

When the self-shielding becomes important, Nn 2 increases and Flw decreases rapidly with r. Then r( c ° o1 ) 
is not much larger than r s h- In such a case, we put r^ co ° ' lower bound. Then r (cool) is given by 

the lesser one of those given by equations ( p5f ) or 

In all the same way as above, we can obtain the value of r( cq ) ; 

r (e q ) =min[r (eq) sh)rsh]) (29) 

where 

(cq) _ r L hW ,1/2 ( ~ Q] 

r no - s h ~ [ (eq) J (M) 
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M M) = ^ nmpr (inf)- (32) 



- 2 0xl0 22 f Llw ) 1/2 ( toB ) 1/2 (31) 

The baryonic mass within the region of influence M^ mf ^ is then 

T 

1.4 X 10 14 M Q (# T )-3/»( 10Mm ^ r _ i& _ O 3/2 (T r g) -6 (ls g =T) -5/4 ( ^ ) -3/4 

0.85 x 10 n M Q ( 10 . 4crg ^ lHz „0 3/2 (3^) 3/2 (T^) }> (33) 

1.0 x lO 8 M (^)- 1 ( loa4CTg ^w iHz _ i )(_T_ 

In equation (p3|), each expression corresponds to the case of r^ inf ^ = r^ ^, r^_ sh , and r s h from the top to 
bottom, respectively. We shall keep this order hereafter. At first glance, the first expression in equation ( |33| ) 
seems to be always larger than the others, but its stronger dependence on temperature makes it important 
for higher temperature (i.e., more massive) objects than the normalized value. From equation (|33|), we can 
see that the mass within a region of influence of an O star already exceeds the scale of the small pregalactic 
object. 

We have considered in this letter the regulation of star formation by photodissociation of molecular 
hydrogen in a pregalactic cloud. On the other hand, Lin & Murray (1992) considered only the regulation 
by photoionization. In such a case, the mass affected by an OB star, namely the baryonic mass within a 
Stromgren sphere, is 

r (St) _ m p nQ 



M w = (34) 



where Q* is the flux of ionizing photons by a OB star and Q* ~ 10 49 s _1 for an 05 star. This is by far 
smaller than our estimated mass of photodissociative influence M^ lnf ' . 

To be more specific in the cosmological context, we consider here pregalactic clouds at virialization. 
The number density at virialization is 

n vir = 0.68cm- 3 ^ (^)(i±i^) 3 , (36) 

and the virial temperature is 



6.8 x l0 2 K^ 3 (^)- 2 / 3 (^^) 2 / 3 (i±^). (37) 
y 0.Q5 J y 10 4 M Q ' v 30 



Substituting equations (36) and ( |37| ) into equation (|33[), we obtain 

2.3 x lO 15 M 0(#T )- 3 / 2 ( lo .^ 

Mf f) = min { 0.58 x 10"M Q ( 10Mm *™. 1H ,. ^ (^^lo^H^T ) 3 

2.2 x 10aM Q(#T )-H 10 . 4crg ^H Z -^ 

(38) 

In order for the star formation to continue after a massive star forms, the region of influence must be 
smaller than the original pregalactic cloud. Then M^ nl ^ < Mb is the necessary condition, which leads to 



1.9 x 10 6 M Q(# ,r 3 / % ^ r ^ ) 
0.58 x 10 11 M Q ( 10 ^f_ lHz _O 3 / 2 (^) 3 / 2 fe 2 o (^)(l^) 3 
4.0 x lO6M (^)- 3 /5( lo2 , cr ^ 



(39) 
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On the other hand, the baryonic mass of a pregalactic cloud that has virial temperature T v - n is 

■ ^vir n3/2i-1/ ,, 1 + Z 

1000K ; 501 0.05 n 30 



M h = 1.8 x lO 4 M (^)3/^-(^)(i±^)-3/ 2 . (40) 



Comparing equations ( p9[ ) and ([10|), we can see that for a small pregalactic cloud {T v -„ < 10 K) the 
condition M^ a1 ' < Mb is hardly satisfied in the redshift range of 10 <C z V j r < 100. This indicates that 
FUV radiation from one or a few OB stars prohibits the whole small pregalactic cloud from H2 cooling and 
quenches subsequent star formation in it. 

After the death of the first OB star, star formation could occur somewhere in the cloud, and another 
OB star could form successively. Thereafter, some massive stars might form one after another, but only 
a few could co-exist simultaneously, as we have shown above. The timescale for reformation of OB stars 
depends on that for H2 replenishment after the death of the dissociating OB star, which is 

Wcool) 

*rep = { (41) 

= 3xl °V<T^)- a ' 2 <Ii?K)- 4 (T^)-'<^|)-" 3 . («) 

If the ionization degree is as high as unity, which is typical in the HII region, this timescale can be very 
short, namely, t rcp ~ 3 x 10 4 yr. If the timescale for reformation of OB stars is smaller than the lifetime 
of OB star, it would be possible for a few OB stars to form every few million years, and, as a result, 
a considerable amount of stars would form in a Hubble time. However, we do not expect that the star 
formation continues as long as a Hubble time, because the gravitational binding energy of baryonic gas in 
such a small pregalactic cloud, 

E„ = °-f^ (43) 

iXvir 

~ 3.3 x 10 48 ergs(^^) 5/3 (— )" 2/3 ^ 3 (^ ±i ^), (44) 
s 4q 4 M ; V 0.05 ; 50 v 30 ; ' v ' 

where M is the total mass, and i? v ; r is the virial radius of the cloud, is so small that a few supernova 
explosions would blow out such a small pregalactic cloud (see, e.g., Dekel & Silk 1986). Thus, our probable 
scenario is as follows: in a small pregalactic cloud, star formation occurs only in a photodissociatively 
regulated fashion until several supernovae explode, and it stops thereafter. 

If numerous small stars were formed per a massive star, a substantial proportion of the original cloud 
would be converted to stars at last. However, in the primordial circumstance, the stellar initial mass 
function could be strongly biased toward the formation of massive stars because of the higher temperatures 
relative to the present-day counterpart. If this was the case and OB stars were formed selectively, the 
amount of gas mass that was converted to stars would be extremely small. 



3. Summary 

We have studied the H 2 photodissociation region around an OB star in a primordial gas cloud. A 
region as large as the whole small pregalactic cloud is affected by only one or a few OB stars and becomes 
unable to cool in a free-fall time under the condition appropriate to virialization. Therefore, in those 
clouds which have virial temperatures less than 10 4 K, star formation does not occur efficiently, unless the 
primordial initial mass function is extremely weighted toward low mass stars. 
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If the rcionization of the universe is caused by stellar UV radiation, some OB stars must form. However 
as we have shown, an OB star formed in a small pregalactic cloud would inevitably photodissociate the 
whole cloud and subsequent star formation would be strongly suppressed. Therefore, stellar UV radiation 
from small pregalactic clouds cannot play a significant role in the reionization of the universe. 



The authors would like to thank Toru Tsuribe for fruitful discussions, Humitaka Sato and Naoshi 
Sugiyama for continuous encouragement, Evan Scannapieco for checking the English, and the referee, Zoltan 
Haiman, for a careful reading of the manuscript and for useful comments. This work is supported in part 
by the Grant-in-Aid for Scientific Research on Priority Areas (No. 10147105) (R.N.), and the Grant-in-Aid 
for Scientific Research from the Ministry of Education, Science, Sports and Culture, No. 08740170 (R.N.). 



- 9 - 



REFERENCES 

Ciardi, B., Ferrara, A., & Abel, T. 1998, ApJ, 516 

Couchman, H. M. P., & Rees, M. J. 1986, MNRAS, 221, 53 

dc Jong, T., 1972, A&A, 20, 263 

Dekel, A., & Silk, J. 1986, ApJ, 303, 39 

Drainc, B. T., & Bertoldi, F. 1996, ApJ, 468, 269 

Fukugita, M., & Kawasaki, M. 1994, MNRAS, 269, 563 

Haiman, Z., Thoul, A. A., & Locb, A. 1996, ApJ, 464, 523 

Haiman, Z., Rees, M. J., & Loeb, A. 1997, ApJ, 476, 458; erratum 484, 985 

Haiman, Z., & Locb, A. 1997, ApJ, 483, 21 

Hutchins, J. B. 1976, ApJ, 205, 103 

Kepner, J., Babul, A., & Spergcl, D. 1997, ApJ, 487, 61 

Lin, D. N. C, & Murray, S. D. 1992, ApJ, 394, 523 

Martin, P. G., Schwarz, D. H., & Mandy, M. E. 1997, ApJ, 461, 265 

Gnedin, N. Y., & Ostriker, J. P., 1997, ApJ, 486, 581 

Tegmark, M., Silk, J., Rees, M. J., Blanchard, A., Abel, T., & Palla, F. 1997, ApJ, 474, 1 



This preprint was prepared with the AAS IAT^X macros v4.0. 



